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ABSTRACT 

This paper presents a method to reconstruct the position of the line of sight and ori- 
entation of the focal plane of the Planck satellite, which may be expressed by the pa- 
rameters of the effective boresight and roll angles, respectively. The accuracy to which 
these parameters must be determined to avoid compromising the reconstruction of the 
Ci power spectrum is assessed. It is expected that the pointing reconstruction will be 
performed using point source transits in the highest frequency channels, given their 
smaller beam sizes. Estimates are made for the number of galactic and extragalactic 
point sources visible by Planck in these channels. The ability of the method presented 
here to reconstruct the pointing parameters to the required accuracies using these 
point sources is investigated. It is found that the scanning strategy employed and the 
actual orientation of the focal plane both independently influence the available num- 
ber of suitable point sources, from which this method may successfully evaluate the 
pointing parameters. While this paper focuses on the pointing reconstruction of the 
Planck satellite, elements of the analysis presented here may be of use elsewhere. 
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1 INTRODUCTION 

Planck is a European Space Agency satellite designed to 
produce high-resolution temperature and polarisation maps 
of the CMB. It possesses detectors sensitive to a wide range 
of frequencies from 30 to 857 GHz, split between two in- 
struments the HFI and LFI, the high and low frequency 
instruments, respectively. Planck is scheduled for launch in 
2007, when it will be inserted into a Lissajous orbit around 
the second Lagrange point of the Earth-Sun system. Planck 
will spin about its axis once per minute, sweeping its de- 
tectors along an almost great circle, the line of sight being 
almost perpendicular to the spin axis. The spin axis will 
nominally be repositioned every hour, and the roughly 60 
or so circles corresponding to a single spin axis positioning 
may be binned together to form a ring. The spin axis passes 
through the centre of the solar panels and is directed to- 
wards the sun, thus maintaining the rest of the satellite in a 
cone of shadow produced by the solar panels. The scanning 
strategy is determined by the frequency and the positions 
of the spin axis repointings over the course of the mission. 
The baseline scanning strategy is such that the spin axis 
remains in the ecliptic plane throughout the course of the 
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mission and it maintains its sun-wards direction by hourly 
respositionings of 2.5'. In reality this scanning strategy will 
not be employed due to the lack of coverage at the ecliptic 
poles; the actual scanning strategy, however, has yet to be 
finalised. 

This paper will investigate some of the issues concerning 
the evaluation of parameters required for the pointing recon- 
struction of the Planck satellite, specifically those resulting 
from non-zero off-diagonal elements in the inertia tensor as 
defined in the nominal satellite reference frame. The slight 
misalignment between the nominal spin axis of the satellite 
and the nearest principal axis of the inertia tensor defines 
the effective boresight angle and the rotation of the focal 
plane. These parameters may only be established using the 
science data, and they will vary over the course of the mis- 
sion due to changes in the inertia tensor as consumables are 
depleted. It is therefore important to assess what the accu- 
racy requirements on these parameters are, and whether it 
will be possible to extract the parameters to those required 
accuracies on a sufficiently regular basis. 

The variation of the misalignment over the course of the 
mission is not yet well defined and is currently under investi- 
gation. It is, however, expected that the frequency at which 
the roll and effective boresight angles must b e evaluated will 
be of the order of once a week 
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Figure 1. The figure shows the unsubtractable errors on the C<s 
which are dominated by cosmic variance to the left and the finite 
beam size to the right. The solid grey curves enclose the region 
corresponding to the error on each individual multipole, whereas 
the dashed grey curves correspond to the error on multipole bins 
of width 50. The black curves correspond to the additional errors 
on the reconstruction of the C^s due to unknown random pointing 
errors. A beam of FWHM 5', with a noise of 10.3/^K per beam 
and a fractional sky coverage of 0.7 have been assumed. 

In Section |5] the accuracy requirements of the effec- 
tive boresight angle and the roll angle, the rotation of the 
focal plane, are discussed. Methods of evaluating the roll 
and effective boresight angles through the use of galactic 
and extragalactic point sources are discussed in Section |^1 
The number of both galactic and extragalactic point sources 
likely to be detectable by Planck axe discussed in Section 0] 
The simulated data on which our methods are tested is dis- 
cussed in Section |^| the results of which are discussed in 
Section HJ 



2 ACCURACY REQUIREMENTS 

The effect of unknown random pointing reconstruction er- 
rors on the recovered Ce values may be seen in Figure^ This 
figure has been replicated from one first shown bv | Lammarel 
feOOll) , using equation modified from iKnoxl Jl995f) : 



C e V (21 + 1) x f shy V + CiW e ) (1) 

where uj is the noise per beam, f s k y is the fractional sky cov- 
erage and We is the window function. The sensitivity of an 
experiment to the £ multipoles is given by its window func- 
tion. For a map made with a Gaussian beam, the window 

_p2 2 

function is given by: We = e b . If the beam sigma, Ob, is 
known, the decreased sensitivity to the higher £ multipoles 
may be accounted for by multiplying by the window func- 
tion. It is this which results in the exponential increase in 
the noise for the high multipoles in Figure as the noise is 
undiminished by the beam size. For low multipoles the noise 
is dominated by the first term in equation^which expresses 
the sample and cosmic variance. 



Table 1. The requirement on the random pointing error versus 
frequency channel, given fractional sky coverage, f a ky 
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100 


9.2 
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0.23 
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113 


7.1 
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0.18 


0.16 


217 
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10.3 


0.16 


0.15 


353 


5.0 


40.9 


0.29 


0.26 


515 


5.0 


46.3 


0.31 


0.28 



Figure^uses the expected beam size and noise levels for 
the 217 GHz frequency channel and assumes 70 per cent sky 
coverage. The grey curves enclose regions corresponding to 
the unsubtractable noise, given by equation^ The solid grey 
curves correspond to the error on each individual multipole, 
while the dashed grey curves the error on multipole bins of 
width 50. 

The black curves, in Figure correspond to the ad- 
ditional errors in the reconstruction of the C{S due to the 
effects of unknown random pointing errors. This results in 
a decreased sensitivity to the higher £ multipoles as is seen 
in Figure Q The pointing uncertainty results in an effec- 
tive smearing of the beam and hence an increase in the ef- 
fective beam width. If this is well known, a new window 
function may be evaluated and the effects of the pointing 
errors will not be as dramatic as those shown. However, in 
reality the pointing errors will be correlated, smearing the 
beam more in one direction, which will lead to asymmetric 
effective beams. The effect of smearing the beam predomi- 
nately in one direction is to increase the sensitivity to higher 
I multipoles as compared to the uniform smearing shown in- 
Figure^ In other words, for a given magnitude, correlated 
pointing errors will lead to a smaller reduction in sensitivity 
at high I, than random pointing errors. However, given that 
the extent of the correlation of the pointing errors is not 
known and that the reduction in the maximum allowable 
pointing error is only significant for highly correlated point- 
ing errors, the pointing errors are assummed to be random 
in the subsequent analysis. 

Ideally, the pointing accuracy should be such that the 
additional errors introduced are less than the unsubtractable 
noise. 

For each frequency channel, with the exception of 857 
GHz in which the Ces are not observable, the unsubtractable 
errors on the Ces, given the expected beam size and noise 
levels, were evaluated. Hence, the ideal pointing accuracy 
required for each frequency channel may be found. These 
are shown in Table for the cases of 70 and 100 percent 
sky coverage. 

The uncertainty in the roll angle should be such that 
the corresponding uncertainty in the position of the detec- 
tor is less than the requirement set by the analysis shown 
in Table Q The required accuracy of the roll angle is deter- 
mined by a combination of the maximum allowable pointing 
error and the location of the detectors in the focal plane. 
The further from the centre of the focal plane the more the 
detectors will be affected by the roll angle. The required 
accuracy of the roll angle for each frequency channel is dis- 
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Table 2. The requirement on the error in the roll angle, given 
fractional sky coverage, f s ky 



Frequency Channel 
(GHz) 


Roll Anf 

fsky = 0.7 

(') 


;le Error 
fsky = 1-0 
(') 


100 


6.9 


6.3 


143 


4.9 


4.4 


217 


(5.7 


6.2 


353 


8.1 


7.2 


545 


8.3 


7.5 



played in Table |21 again for both the cases of 70 and 100 
percent sky coverage. 

As shown in Table [5] the 143 GHz channel sets the 
requirement on the uncertainty in the roll angle to be 4.4'; 
assuming no other pointing errors from other sources. 



3 METHODS 

In this section methods for evaluating the roll angle and the 
effective boresight, using measurements of point sources, will 
be considered. The effect of a roll angle will be to change the 
positions of the detectors with respect to the scan and cross- 
scan directions, so any method to evaluate the roll angle 
will depend on the focal plane layout. Figure [5] illustrates 
the layout of the focal plane. The four highest frequency 
channels have offset detectors to ensure full sampling in the 
cross-scan direction. The centres of adjacent detectors are 
offset by 1.5' in the nominal cross-scan, with every detector 
being in a pair with the same nominal cross-scan location. 
The highest frequency channels will be the most suited to 
extracting the roll angle and effective boresight angle, due 
to the smaller beam size and the larger numbers of point 
sources which will be detected at the higher frequencies. For 
this paper, unless otherwise stated, we will consider only the 
857 GHz detectors. We will also define a positive roll angle 
as a clockwise rotation, and a negative rotation as any anti- 
clockwise rotation about a fiducial reference point, FRP, in 
the centre of the focal plane. 

In the nominal case of zero roll angle there are two 
pairs of detectors with the same cross-scan position. Any 
separation found between the cross-scan positions of these 
detectors may be used to infer the roll angle, as shown in 
Figure 01 

In order to evaluate the roll angle a method is needed 
to find the cross-scan positions of the detectors, or the cross- 
scan separation of the detectors aligned in the nominal scan 
direction, directly from the amplitudes of a point source as 
observed by each detector. Once a method has been found it 
may be investigated as to what signal-to-noise ratio, SNR, is 
required to obtain the roll angle to the required accuracy and 
hence an estimate of the number of available point sources 
may be made. Methods for reconstructing the roll angle are 
discussed in Section RllI 

The effective boresight is the angle between the spin 
axis and the scan circle which passes through the FRP, 
shown as a in Figure HJ In order to evaluate the effective 
boresight the cross scan position of a known point source 
relative to the FRP must be found. Once the roll angle has 





2 






.143 


.143 


.143 




.143 












.143 


.143 


.143 




.143 




CO 


1 


















CD 
0) 






545 


545 fi57 


£57 


£57 


£57 


545 


545 ; 


(deg 




















c 
o 

o 
</) 







353 


.353 353 


.353 + 
217 


353 

217 


353 


353 
217 


353 _ 


o 




















U 






















-1 






217 
.100 

, r . 


217 

.100 




217 
.100 


2 


7 

.100 





-2-10 1 2 

Scon (degrees) 



Figure 2. This figure displays the focal plane layout of the HFI, 
in terms of the nominal scan and cross-scan directions. The posi- 
tive scan direction indicates the direction in which the focal plane 
traverses the sky. A source will therefore cross the focal plane from 
right to left in this figure. The exact location of the detectors are 
given by the small filled circles. The top four frequency channels 
have offset detectors to ensure full sampling in the cross-scan di- 
rection. The centres of adjacent detectors are offset by 1.5' in the 
cross-scan direction, with every detector being in a pair with the 
same cross-scan location. The roll angle is defined as a rotation 
about the FRP, which is shown by the cross, and the effective 
boresight is defined as the angle between the spin axis and the 
centre of the focal plane. In this paper we will define a positive 
roll angle as a clockwise rotation, and a negative rotation as any 
anti-clockwise rotation about the FRP. 

been evaluated, and hence the cross scan positions of the de- 
tectors are known, the evaluation of the cross scan positions 
of point sources becomes possible. Methods for evaluating 
the cross scan positions of point sources will be discussed in 
Section 13.21 

These methods are applied to data from one or two 
pointings of the spin axis, as required. The data from each 
spin axis pointing may be used in terms of the individual 
revolutions, circles, or the accumulated revolutions corre- 
sponding to the spin axis pointing, a ring. In the text below 
a scan circle may be thought of as the path of the line of 
sight of the FRP on the sky in the case of a circle or the 
average path in the case of a ring. This differentiation is nec- 
essary as in reality the spin axis describes a small rotating 
ellipse on the sky. This nutation of the spin axis results in 
deviations from the expected path of the line of sight given 
a spin axis rotating around a fixed point. These deviations 
are small and will average out in the ring data as the period 
of the nutation is greater than the rotation rate of the satel- 
lite. Providing that the amplitude of the nutation remains 
small it may be safely ignored in all cases, should this no 
longer be the case the effects of the nutation may be cor- 
rected for using information on the nutation parameters as 
found from the Star Tracker data. Throughout this paper 
it is assumed that these methods will be applied to rings, 
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Figure 3. This figure shows the effect of a change in the orien- 
tation of the focal plane with respect to the scan direction. The 
dotted grey circles correspond to detector positions in the nomi- 
nal orientation of the focal plane, whereas the solid black circles 
correspond to their positions after a rotation, p, about the FRP. 
For the highest frequency channels the centres of alternate de- 
tectors are aligned in the nominal scan direction as shown in the 
figure by the dotted grey line. Any difference in the cross-scan po- 
sition of these detectors is therefore indicative of a non-zero focal 
plane rotation. The nominal scan separation, NS sop , of alternate 
detectors may be used with their cross-scan separation, C scp , to 
calculate the roll angle. 

given their lower noise levels and their greater robustness to 
the effects of nutation. 



3.1 Reconstructing the Roll Angle 

3.1.1 Using a single ring 

In the first instance one may consider evaluating the roll an- 
gle using a single ring, using the observed amplitudes from 
the detectors to solve for the amplitude of the point source 
and its cross-scan position. Here a ring consists of the data 
corresponding to a single repointing of the spin axis, and 
the observed amplitude and location in the scan direction of 
the point source is found fo r each detector from either bin- 
ning or by phase-ordering llvan Leeuwen et. al.ll200lh this 
data. Once the amplitude and cross-scan position of the 
point source has been evaluated the amplitude observed by 
each detector may be used to evaluate its cross-scan sepa- 
ration from the point source position. Given these relative 
cross-scan locations the roll angle may be evaluated. How- 
ever, in practice the errors in the detected amplitudes re- 
sult in errors in the cross-scan locations which may lead to 
a large uncertainty in the recovered roll angle. An exam- 
ple of this is shown in Figure [IJ in which the errors in the 
detected amplitude have been converted into errors in the 
cross-scan position of the detector. In this figure it has been 
assumed that the amplitude and cross-scan position of the 
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Figure 4. This figure displays the amplitudes, as observed by 
each detector, for two different roll angles, given a point source 
with an amplitude of 100 times the ltr noise and cross-scan po- 
sition from the FRP of 0?552. The circles represent the detected 
amplitudes given a roll angle of 0?4, while the triangles 0?7. The 
la errors in the amplitudes have been converted into errors in 
the cross-scan position of the detectors. The corresponding ltr 
uncertainties in the positions of the detectors overlap, it is there- 
fore not possible to establish the true roll angle without a higher 
signal-to-noise point source. 



point source have been successfully recovered with negligi- 
ble error; in this example the amplitude of the point source 
is 100 times the la noise and its cross-scan position from 
the scan circle which contains the FRP is 0?552. For com- 
parison two scenarios are plotted; the circles represent the 
detected amplitudes in the four detectors resulting from a 
0?4 roll angle, while the triangles correspond to the ampli- 
tudes resulting from a 0.7 roll angle. The la errors in the 
amplitudes have been converted into errors in the cross-scan 
position of the detectors. Working across the figure pairing 
up the adjacent triangles and circles which correspond to the 
same detector, it is possible to see that uncertainties in the 
cross-scan positions, for each of the roll angles, overlap for 
all of the detectors. Hence in this example it is not possible 
to successfully distinguish between roll angles of 0?4 to 0?7. 
Other cross-scan positions of the point source with respect 
to the FRP result in similar difficulties over different ranges 
of roll angle. 

In reality there will be errors in the recovered amplitude 
and cross-scan position of the point source, together with er- 
rors in the relative calibration between detectors which will 
increase the effective errors in the cross-scan positions of the 
detectors and hence the uncertainty in the recovered value 
of the roll angle. Indeed, it is not possible to evaluate the 
cross-scan position of the point source, without first know- 
ing the value of the roll angle due to the lack of knowledge 
on the positions of the detectors. 
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Figure 5. This figure shows two neighbouring, although not nec- 
essarily adjacent scans, which pass on either side of a bright point 
source. The detector positions relative to the scan are denoted by 
Yj, while the Xij represents the offset from the point source of 
detector j in scan i. d is the cross-scan position of the point 
source from scan i and S12 is the distance between the two scans 
at the location (longitude) of the point source. 



3.1.2 Using Two Rings 

If there are no abrupt changes in the inertia tensor during 
the mission, such as fuel shifts, the roll and effective bore- 
sight angles will only change slowly over the duration of the 
mission. Assuming that this is indeed the case and hence 
the roll angle is a slowly varying parameter, its evaluation 
may be improved by including the additional information 
available from the detection of the same point source on a 
neighbouring ring. Figure shows schematically the rela- 
tive locations of the detectors to the scan directions, cor- 
responding to the two rings, and the point source position. 
The cross-scan position of the point source from scan 1 and 
scan 2 is given by Ci and C2, respectively; C\ has a negative 
sign so that the sign convention for the cross-scan position 
may be the same for both scans. As the roll angle is un- 
changed, the cross-scan positions of the detectors, Yj, from 
the FRP, will therefore not have changed between the two 
neighbouring scans. The sum of the offset from the point 
source position of a detector in scanl, xij, and the same 
detector in scan2, x 2 j will therefore be equal to the sum of 
the angles Ci and C2. The angular separation between the 
two scans at the location of the point source, S12, may also 
be related to the sum of the angles Ci and C2, however in 
this case this relationship is only an approximate. Though 
as will be shown below for neighbouring scans, it may be 
thought of as being exact: 



The angle between the two scans, S12, at the location of the 
point source depends on the separation of the mean spin axis 
positions corresponding to each of the scans, and the eleva- 
tion, ip, of the point source from the plane that connects the 
mean spin axis positions. Figure |S| shows a schematic of the 
two scans and the point source in question. The separation 
between the mean spin axis positions is given by <p, this is 
also the maximum separation between the two scans. The 
angles, ip and 4>, may be evaluated using additional infor- 
mation provided by the Star Trackers. In the case of <f>, it 
may be evaluated from the two recovered mean spin axis 
positions. The evaluation of ip also requires the time of first 
transit of the FRP though the reference circle which con- 
nects the mean spin ax is position and the North Ecliptic 
Pole, NEP, as defined in Ivan Leeuwen et. alJ <l200ll) . 

The angle 5*12, also has a dependence on the effective 
boresight, a, though given the expected value, a « 85°, this 
dependence is weak: 



cos(5i2) = cos{</>) (cos 2 (a) + sin 2 (a) cos 2 (^>)) 
+ sin 2 (a) sin 2 {ip). 



(3) 



As can be seen from Figure |S] using the value of £12 for 
the sum of the two cross-scan angles to the point source, 
— Ci + C2, will result in an underestimation of this value. 
This underestimation will be the greatest when the point 
source lies midway between the two scans. To quantify the 
magnitude of this underestimation, the angle between the 
two cross-scan directions, /3, given a point source which lies 
midway between the two scans may be found using: 



cos(/3) = cos(-7r — </>) (cos 2 {ip) cos 2 (a) + sin 2 (a)) 



cos 2 (a) sin 2 (i/j). 



(4) 



Given that we are only interested in neighbouring scans the 
angle, <f> will be small and the angle, 5, between Ci,2 and 
5*i2 may then be found using: 



S = 



(•») 



For 0^1°, which corresponds to ^ 24 nominal repointings 
of 2.5', the underestimation in the sum of the cross-scan 
angles is of the order of a thousandth of a percent. This 
corresponds to an overestimation for the magnitudes of the 
roll angle considered here of the order of a hundredth of an 
arc-second. Hence, for our present purposes the relation in 
equation |5] may be considered to be exact. 

From Figure it is seen that the sum of the cross scan 
angles C\ and C2 is equal to the sum of the offsets, from 
the point source position, of the j th detector in scans 1 and 
2, xij and X2j- The angle S12 is therefore also equal to the 
sum of these offsets, via equation [5] 



5*12 = Xij + X2j. 



(6) 



The cross-scan position of the point source from the scan 
and the offsets of the detectors from the point source, may 
be related to the position of the detectors relative to the 



— C\=x\ 
C 2 = x 2 j+Yj. 



(7) 
(8) 
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Figure 6. This figure shows the angles discussed in the text relat- 
ing to two neighbouring scan circles with a point source located 
between them. The top diagram shows the plane in which the mo- 
tion of the mean spin axis has occurred, and the angle between 
this plane and the point source, ip. The angle between the spin 
axis and the scan circle, a, is the effective boresight. The angle tp, 
may be evaluated once the reference phase, or time of first transi t 
of the FRP is known, as described in lvan Leeuwen et. alJ 1200 ll) . 
In the first instance this may be evaluated by the Star Trackers. 
The lower diagram shows the angle, <f>, through which the spin 
axis has been repositioned. The angles C\ and C2 represent the 
cross-scan angular distance between the point source and scan 
circles 1 and 2, respectively. The angle S12 is the angular separa- 
tion between the intersection points of C\ and C2 on scan circles 
1 and 2. The angle rf>, between the two spin axis positions, has 
been exaggerated for clarity in this figure; in the analyis presented 
here d> never exceeds f°. 



The observations of the point source may be described by 
the following, the amplitude of the point source as seen by 
the j th detector in the i th scan is denoted by Aij\ 



Ai 



A exp 



2a? 



(9) 



Where en, is the lcr width of the assumed Gaussian beam 
for the j th detector. 

The expression for the ratio of the amplitudes of the 
point source as observed in the two scans by the j th detector 
may be rearranged to find: 



2 

x 2j 



2 

X Xj 



2a j In 



Ail 
A 2j 



from which, using equation |S| we find: 



x 2 j - XX j 



2<7| 

S12 



A 2j 



(10) 



(11) 



where everything on the RHS of this equation is known 
observationally. Equation 1111 shows a major benefit of this 
method; as it is based on the ratio of the amplitudes seen 
by the same detector in different scans, uncertainties in 
the relative calibrations between detectors will be irrelevant. 

An alternative expression for the difference between the 
offsets of the j th detector from the point source on each scan, 
may be found using equations [7| and [S] 
1 



((C 2 + Ci) - (x 2j - . 



(12) 



The cross-scan separation, C scp , of the aligned detectors 
may be found: 

1 
2 
1 



y 3 -Yi 



((221 - x-n) - (£23 - £13)) • 



Y 4 -Y 2 = ^ ((x 2 2 - X12) - {x 2 4 - X 14 )) 



C s . 



Y 3 -Yi_ = Ya-Y 2 , 



(13) 

(14) 
(15) 



and hence the roll angle, p, may be evaluated, where the 
denominator, NS sep , is determined by the separation in the 
scan direction of the aligned detectors in the case of zero roll 
angle, and is shown in Figure [!|] 



3.2 Evaluating the Effective Boresight 

The effective boresight, which is the angle between the mean 
spin axis position and the scan circle, may be found if a point 
source of known position is observed in the data from the 
scan circle. The cross scan position, C, of the point source 
is found with respect to the FRP, and hence the scan circle. 
Here the cross scan position of the point source is defined as 
positive if the point source lies on the opposite side of the 
scan circle to the spin axis position. The effective boresight, 
a, may then be found using: 



C, 



(17) 



where uo is the angle between the point source position, 
(f3 p t,\pt), and the mean spin axis position, (f3 sa ,\ S a), and 
is given by: 
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cosuj = cos f3 pt cos/3 sa cos(A pt - X sa ) + sin f3 pt sin f3 sa ■ (18) 

The cross scan position of a point source may be found 
using data from a single ring or by an extension of the two 
ring method for the evaluation of the roll angle. 

In the case of data from a single ring, the amplitudes 
found from the transits of the point source are fit by a Gaus- 
sian beam profile to establish best fit values for the position 
and amplitude of the peak. This method will only provide 
unbiased cross scan position for the point source if it lies 
between the range of cross scan positions of the detectors. 
Hence the value of the roll angle is required prior to the 
evaluation of the cross-scan position of the point source. 
Additionally, uncertainties in the relative calibration and/or 
beams of the detectors will increase the uncertainty in the 
recovered value of the cross-scan position and hence effective 
boresight angle. 

Alternatively, by an extension to the method to extract 
the roll angle from two rings, the cross-scan position of the 
point source may be found using the ratio of the amplitudes 
observed by each detector in two neighbouring scans. 

The sum over all detectors using equation 1121 may be 
used to evaluate the difference between the cross-scan posi- 
tions of the point source to each scan: 

N N 

3 i 

where N=4 using a single frequency channel, or N=8 com- 
bining both the 857 and 545 GHz channels. 

Once the roll angle is found the positions of the detec- 
tors, Yj, are known, indeed the sum over the positions of the 
detectors is approximately a constant for small values of the 
roll angle. The cross scan position of the point source may 
be found for each scan by combining eg uations 121 and 1 1 91 

a N / / \ \ N 

Once the cross scan position of the point source is 
known with respect to the FRP, the effective boresight may 
be evaluated via equations 1171 and 1181 if the position of the 
point source is known. 

4 PREDICTING THE NUMBER OF POINT 
SOURCES VISIBLE WITH PLANCK HFI 

The numbers of point sources visible with Planck HFI may 
be pr edicted by using the I RAS point source catalogue 
(PSC, Bcichma n et. alJ (|l98S|)). This is a catalogue of some 
250,000 well-confirmed point sources, providing positions, 
flux densities at 12, 25, 60 and lOO^m, uncertainties and var- 
ious cautionary flags which are given for each source. Only 
the brightest of these sources will be visible by Planck. Since 
the brightest of these sources have an increased chance of 
actually being extended objects, it is important to establish 
the characteristics which best indicate a point source na- 
ture. In order to establish these characteristics a population 
of objects which are known to be point sources is required. 



Such a population of point sources may be found by apply- 
ing th e following stellar colour cuts llRowan- Robinson et. alJ 
1986) to the PSC: 

log(ff) < 0-2- (22) 

The resultant population of stellar, and hence point 
sources may then be investigated. The population consists of 
450 bright, S\oo ^ 2.0 Jy, point sources which were all found 
to have a 60u-m correlation coefficient of grade A, where the 
correlation coefficient is an indication of how closely the data 
match the template for response from a point source. The 
100|J.m correlation coefficient was found to be grade A for 
331 of these sources, grade B for 104, and grade C for 12. 
The reliable indicators of a point source were then taken to 
be a 60 (xm correlation coefficient of grade A and a 100u.m 
correlation coefficient of grade A or B. 

The above cuts on the correlation coefficients were ap- 
plied to the IRAS PSC, excluding those sources with correla- 
tion coefficients worse than grade A at 60u.m and grade B at 
lOOu-m. In addition sources with only upper limits at 100 or 
60u.m were also excluded. The colour cut log (4^) > —0.3 
was applied, in order to remove all the sources with rapidly 
declining fluxes towards longer wavelengths. After these cuts 
the number of sources remaining is 7874, 558 8 of which are 
in the PSCz catalogue iSaunders et. al.ll2000Tl . the rest pre- 
dominately populate the galactic plane. 

The IRAS galax ies may be modelled using the work of 
iDunne et. alJ ( 2000) and IDunne fc Ealesl (1200 ll) who made 
850u.m and 450 u.m maps of galaxies fr om the IRAS Bright 
Galaxy Sample. IDunne fc Ealesl (J2001) find that this com- 
bined data is best fit by a two component spectral energy 
distribution, SED. The two components are modelled as two 
/3 = 2 grey bodies corresponding to warm and cold dust tem- 
peratures: 

Su = N w v p B(v,T w ) + N c u B(v,T c ), (23) 

where B(^,T) is the Planck function for both components, 
T m and T c are the temperatures of the warm and cold com- 
ponents and N w and N c correspond to the relative masses 
of the dust in each component. 

The extragalactic point source fluxes may now be ex- 
trapolated to Planck frequencies. Since a two component 
dust model is a four parameter fit, two temperatures and two 
normalisations, and as the 12(im flux is only present as an 
upper limit for so me of these sourc e s, the average cold tem- 
perature found by IDunne fc Ealesl (1200 if) of 20. 9K is used 
and the remaining three parameters are solved for using the 
25, 60 and 100u.m fluxes. Using these parameters the fluxes 
of each galaxy at 350, 550 and 850u.m , corresponding to the 
Planck frequencies of 857, 545 and 353 GHz, are found. 

The number of point sources above a required flux level 
at each frequency may now be found. The flux levels cho- 
sen are those which correspond to a range of signal-to-noise 
ratios in a binned ring. This allows a direct comparison 
between frequency channels as each frequency channel is 
expected to have different noise levels. Table [3] shows the 
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Table 3. Number of extragalactic point sources above a flux 
which corresponds to a given SNR in a binned ring 



SNR 






Frequency (GHz) 








857 




545 




353 






S (Jy) 




S (Jy) 




S (Jy) 




10 


2.84 


452 


2.51 


82 


1.54 


26 


50 


14.20 


59 


12.55 


14 


7.70 


4 


100 


28.40 


22 


25.10 


7 


15.40 


1 


200 


56.80 


9 


50.20 


4 


30.80 





400 


113.60 


5 


100.40 





61.60 






Table 4. Number of galactic point sources above a flux which 
corresponds to a given SNR in a binned ring 



SNR 






Frequency (GHz) 








857 




545 




353 






S (Jy) 


N^S 


S (Jy) 


N>S 


S (Jy) 


N^S 


10 


2.84 


1147 


2.51 


575 


1.54 


269 


50 


14.20 


545 


12.55 


203 


7.70 


69 


100 


28.40 


361 


25.10 


122 


15.40 


33 


200 


56.80 


219 


50.20 


57 


30.80 


14 


100 


113.60 


131 


100.40 


25 


61.60 


4 



number of extragalactic point sources above a given signal- 
to-noise ratio with the corresponding flux for the top three 
Planck frequency channels. 

The galactic point sources are modelled using a single 
grey body fit using (5 = 2.3, as expected for compact HII 
regions, M-type stars and Carbon stars, to extrapolate to 
Planck frequencies: 



p p 2hv 



(24) 



c 2 exp(hv/kT) - 1' 

Table [I] shows the number of galactic point sources 
above a given signal-to-noise ratio with the corresponding 
flux for the top three Planck frequency channels. 

As seen in Tables [3] and 2] the vast majority of the point 
sources will be detected in the highest frequency channel of 
857 GHz. Once the number of point sources above given 
signal-to-noise ratios have been found, it is possible to es- 
timate the number of them from which the roll angle may 
be evaluated. Not all potential locations of a point source 
relative to two neighbouring scans, shown by Figure may 
successfully be used to extract the roll angle as will be ex- 
amined below in Section |5] 



5 SIMULATED DATA 

In order to test the performance of this method, simulated 
observations were generated. This data was comprised of the 
peak amplitudes of the transit of a point source for each of 
the four 857 GHz detectors in two neighbouring scans. 

In generating the simulated observations, the location of 
the point source was fixed to lie on or between two successive 
scans, separated by a single repointing, shown schematically 
by the shaded area in Figure |7] The signal-to-noise ratio of 
the point source, its position between the two scans and the 
angle between the point source and the plane of motion of 



the spin axes, ip in Figure |S| were varied for different roll 
angles. Henceforth the angle, ip will be referred to as phase. 
It should be noted that this defin ition of phase i s diffe rent 
from that which appears in Ivan Leeuwen et. alj 11200 II) . In 
the simulated observations it is assumed that the separation 
between successive repointings of the spin axis is the nom- 
inal 2.5'. The peak amplitudes of the transit of the point 
source are then found in two neighbouring scans, separated 
by a given number of repointings. In summary the simulated 
observations vary the following parameters: 

(i) The cross-scan location of the point source between 
the two successive scans 

(ii) The phase, tp, of the point source 

(iii) The signal-to-noise ratio of the point source 

(iv) The number of repointings between the neighbouring 
scans 

(v) The roll angle 

By varying the phase, the signal-to-noise ratio and the 
position of the point source between the two successive 
scans, the dependence of these parameters on the ability of 
this method to evaluate the roll angle may be investigated. 
The results presented here are from 200 noise realisations 
for every combination of the parameters investigated. Un- 
less otherwise stated the error assumed in the mean spin 
axis position was 1" and the error in the reference phase 
was 20'. This error in the reference phase corresponds to 
the worst case initial value from the Star Tracker. 



6 RESULTS 

6.1 Recovering the roll angle 

Two criteria must both be met before this method may be 
used to evaluate the roll angle, 

(i) The point source is detected at or above the 5a limit 
in all detectors in both scans 

(ii) The point source is found to lie above the detectors 
in scan 2 and below the detectors in scan 1, as shown in 
Figure |S] 

These two criteria, together with Tables|3]and2] may be 
used to show that the 545 GHz channel may not be of prac- 
tical use in extracting the roll angle. The 545 GHz detectors 
are further from the FRP and are hence more affected by 
the roll angle than the 857 GHz detectors. For non-zero roll 
angles the separation of the two neighbouring scans, which 
meet the constraint of our method that the point source is 
required to be above all detectors in scan 2 and below all 
detectors in scan 1, rapidly becomes wider than 6a b, where 
at is the la width of the assumed Gaussian beam. This re- 
quires a point source with a flux > 120 Jy in order for the 
first criterion, the detection of the point source at or above 
the 5cr level in all the detectors in both scans, to be met. 
As seen in Table [3] there are only 25 galactic point sources 

> 100 Jy at 545 GHz, and there are in fact no point sources 

> 120 Jy expected at 545 GHz. Given the low numbers of 
point sources available at 353 GHz, it is therefore expected 
that the roll angle will be evaluated with data solely from 
the 857 GHz channel. 

Using data from the 857 GHz channel, the two succes- 
sive scans are considered first and if the two criteria are met, 
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Figure 7. This figure shows a schematic of several neighbouring 
scans. In the simulated data the point source is constrained to lie 
on or between two adjacent scans, corresponding to the shaded 
area in this figure. The roll angle is evaluated using the closest 
possible two scans in which the requirement that the point source 
is located above the detectors in scan 2 and below the detectors in 
scan 1, is met. The point source must also be sufficiently bright to 
be detected at the 5c level in all the detectors on the two scans. 
It can be seen that the ability to evaluate the roll angle cannot 
naively be related to the signal-to-noise ratio of the point source, 
but will also depend on the relative location of the point source 
to the scan circles and on the roll angle itself. 

then an attempt may be made to reconstruct the roll angle. 
If the first criterion is met, but the second is not then the 
next closest scan to the point source is considered and the 
status of the above criteria are reaccessed. The considera- 
tion of further repointings on either side of the shaded area 
in Figure Q continues until both criteria are met, and the 
roll angle may be evaluated, or the first criterion is failed. 

As the signal-to-noise ratio of the point source increases, 
the detections will be above the 5a limit in two scans sepa- 
rated by a larger number of repointings. For scans separated 
by a larger number of repointings the point source is more 
likely to meet the second criterion. In other words, as the 
signal-to-noise ratio increases the number of cross-scan and 
phase locations for the point source, in the shaded region in 
Figure|7| which allow an attempt to be made in reconstruct- 
ing the roll angle, increases. This may be seen in Figure [S] 
which shows the actual area available to this method versus 
the phase of the point source for several different signal-to- 
noise ratios. The shape of the curves in Figure |S| are de- 
termined by the two criteria which must be met. In reality 
both criteria may be met for a number of scans, in this case 
the scans separated by the largest number of pointings will 
correspond to the greatest coverage of the shaded region in 
Figure |7| It is this maximum value for the actual area avail- 
able which is plotted in Figure |H| The dependence of the 
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Figure 9. This figure shows the number of repointings between 
the two scans used versus the phase of the point source for zero 
roll angle and a 100 signal-to-noise point source. The grey curve, 
shown for comparison, is the area available to the method for a 
100 signal-to-noise point source and zero roll angle as previously 
shown in Figure 181 

curves in Figure [£] on the two criteria may be more easily 
shown in Figure [5] In this figure the number of repointings, 
of the nominal 2.5', between the two scan circles used is plot- 
ted against the phase of the point source. The area available 
to the method for a 100 signal-to-noise ratio point source is 
overplotted for comparison. 

As the phase increases, the angular distance between 
the two scan circles decreases, this increases the signal-to- 
noise of the detections of the point source in each scan. It 
is this effect which is the cause of the first increase in the 
area available in Figure [5] which is seen to occur when the 
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Figure 10. This figure shows a contour plot of the area available 
to this method given a point source with a signal-to-noise ratio 
of 100 and zero roll angle. The darker regions correspond to a 
larger available area. The curves in Figures Isl and l9l correspond 
to the maximum value of the available area for a given value of 
the phase. The black/white dots correspond to the value of phase 
and spin axis separations used in the simulated data. 



number of repointings between the two scans used increases. 
This is due to the decreasing angular separation between 
the scans which has allowed the first criterion to be passed 
by two scans separated by an additional repointing. This 
additional repointing increases the separation between the 
detectors on each scan used and thus increases the chance 
of the second criterion being met. The continued increase 
in the available area with a constant number of repointings 
between the scans used is again due to the decreasing an- 
gular separation of the scans which increases the area over 
which the first criterion will be met. The subsequent down- 
turn in the area available is due to the effect of the second 
criterion, as the separation of the detectors on the two scans 
is decreasing there are fewer possible locations for the point 
source which allow it to be positioned above the detectors in 
scan 2 and below them in scan 1. This process of increasing, 
then decreasing area available with phase is repeated once 
the angular separation between the scans decreases enough 
to allow the first criterion to be passed on two scans sepa- 
rated by a larger number of repointings. This cycle of varia- 
tion with phase continues until at large phases the total area 
between two successive scans is available to this method. 

Figure 1101 shows a contour plot of the area available, 
with the darker regions corresponding to a greater available 
area. The curves in Figures|S]and|U]correspond to the max- 
imum value of the available area for a given value of the 
phase. The black and white dots correspond to the value of 
phase and spin axis separations used in the simulated data. 
The shape of the contours are defined by the two require- 
ment criteria; the contours towards narrower spin axis sepa- 
rations are governed by the second criterion and the contours 
towards wider spin axis separations are determined by the 
first criterion. 
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Figure 11. This plot shows the actual area available to the 
method described here versus the phase of the point source, of 
a signal-to-noise ratio of 100, for different roll angles. The grey 
curve shows the actual area enclosed by the two scans in each 2? 
phase bin 



As the roll angle changes the relative positions of the 
detectors to the scan, it will also affect whether the second 
criterion is met. Figure fTTI shows the actual area available 
versus the phase of a point source, of signal-to-noise ratio 
100, for several different roll angles. The differences between 
positive and negative roll angle may be explained by the fact 
that the positions of the detectors in the focal plane are not 
symmetric about the FRP. A negative roll angle increases 
the area covered in the cross-scan direction by the detectors 
more than a positive roll angle of the same magnitude. The 
larger the area covered in the cross-scan direction the smaller 
the chance of the second criterion being met. 

Figure [T2l shows how the percentage area of the shaded 
area in Figure [7| which corresponds to the area between two 
successive scans, varies with the roll angle for point sources 
of different signal-to-noise ratios. The area available to this 
method is determined by the requirement criteria; (i) the 
point source detections in each scan being above the re- 
quired signal-to-noise ratio and (ii) the location of the point 
source with respect to the detectors in each of the scans 
used to extract the roll angle. The grey curves shows how 
the percentage area varies with roll angle using the input to 
the simulation to access whether the position of the point 
source meets the second criterion. The black curves use a 
routine to evaluate whether the point source lies above or 
below the scan in question based on the relative amplitudes 
of the point source as seen by each detector. This routine 
makes no assumptions on the value of the roll angle for this 
evaluation. While this will initially be the case, once the 
roll angle has been evaluated and its time variability mod- 
elled, the routine to evaluate whether the second criterion is 
passed may incorporate the previous values of the roll an- 
gle to improve upon its acceptance/rejection of the use of a 
point source for the roll angle evaluation. 

Figure tT^l shows the mean error on the recovered roll an- 
gle versus the actual value of the roll angle for point sources 
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Figure 12. The plot shows how the percentage of the area avail- 
able varies with the roll angle, for point sources of varying signal- 
to-noise ratios. The black curves correspond to the number of 
point sources which are found to meet the two requirement cri- 
teria, discussed in the text. Whether the second criterion is met 
is assessed using a function which determines whether the point 
source lies above or below the detectors on a scan based on the 
ratio of their observed amplitudes. The grey curves correspond to 
the number of point sources included in the simulated data which 
meet the two requirement criteria. 
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Figure 13. This plot shows the error on the recovered roll angle 
against the actual roll angle for point sources with different signal- 
to-noise ratios. The grey line shows the requirement on the roll 
angle error of 4.4'. 



Figure 14. This plot shows the error on the recovered roll angle 
against the phase location of the point source for varying signal- 
to-noise ratios. The grey line shows the requirement on the roll 
angle error of 4.4'. 




Figure 15. This figure shows contours of constant mean error in 
the recovered roll angle for zero roll angle and a point source with 
a signal-to-noise ratio of 100, against the phase of the point source 
and the angular separation of the spin axes of the scans used. All 
the contours shown are have values below that of the requirement 
on the error in the roll angle, as discussed in Section[3] The darker 
regions correspond to lower recovered errors on the roll angle. The 
black/white point correspond to the actual values of phase and 
spin axis separations used in the simulated data. 



of different signal-to-noise ratios, averaged over the 200 real- 
isations of the noise used and the other parameters listed in 
Section [S] The grey line is the requirement on this error as 
outlined in Section[5] The error on the recovered roll angle is 
not dependent on the actual value of the roll angle for small 
magnitude roll angles, but as the roll angle increases so does 
the error on the recovered roll angle. Indeed, for larger roll 



angles it will be difficult to recover the roll angle using this 
method. This is due not only to the increasing error in the 
recovery of the roll angle, but also to the decreasing area 
available to this method with increasing roll angle as seen 
in Figure H2l This decrease in area available corresponds to 
a decrease in the number of suitable point sources. 

Figure 1141 shows the mean recovered roll angle error, 
over the noise realisations and other parameters investi- 
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gated, for different signal-to-noise ratio point sources versus 
the location in phase of the point source. The increase in 
the value of the mean recovered error at large phases is due 
to the error in the reference phase. This error results in an 
uncertainty in the phase of the point source which results in 
an error in the angular separation, via equation^ of the two 
scans used to recover the roll angle. Given the value used for 
the error in the reference phase is the expected worst case 
scenario for this value, the error in the recovered roll angle 
will have very little dependence on the error in the reference 
phase. 

Figure 1151 shows the dependence of the mean error in 
the roll angle on the phase, and the angular separation of the 
spin axes corresponding to the two scan circles used in the 
evaluation of the roll angle. In this plot a point source with 
a signal-to-noise ratio of 100 was used to reconstruct the 
roll angle, which in this case has the value of zero degrees. 
All the contours in this plot exceed the requirement on the 
error in the recovered roll angle as discussed in Section [5] 
with the darker regions corresponding to lower errors. The 
black/white points correspond to the actual values of phase 
and spin axis separations used in the simulated data. This 
figure shows that optimum phase location for a point source 
is at intermediate to high phase, with the error increasing 
again above ~80° in phase. 

6.2 Recovering the cross-scan position 

Using equations 1201 and 1211 from Section 13.21 it is possible 
to evaluate the cross-scan position of the point source used 
in evaluating the roll angle as above. This method of evalu- 
ating the cross-scan position of the point source is immune 
to relative calibration errors between the detectors and can 
incorporate different beams for each of the detectors. 

Figure 1161 shows the mean error in the cross-scan po- 
sition versus the location in phase of the point source, for 
point sources of different signal-to-noise ratios. The other 
parameters remaining in the simulated data are averaged 
over together with the 200 noise realisations. The increase 
in the error towards larger phases is due to the dependence 
of the cross-scan position on the angular separation, S12, 
between the two scans at the phase location of the point 
source. The error in the angular separation of the two scans 
is dependent on the error in the reference phase; in these 
simulations the uncertainty in the reference phase is taken 
to be the worst case scenario value of 20'. For more realistic 
values of the error in the reference phase the error in the 
cross-scan position at large phases will reduce towards the 
low phase values. 

Figure [T71 shows the mean error in the recovered cross- 
scan position versus the roll angle, for different signal-to- 
noise ratio point sources. This figure shows that the error in 
the recovered cross-scan position follows the same behaviour 
as the error in the recovered roll angle; no dependence on the 
actual value of the roll angle for small magnitudes, with an 
increase in the error towards larger magnitude roll angles. 
Since the other parameters in the simulated data have been 
averaged over, the larger errors in the recovered cross-scan 
position at high phases, as seen in Figure Hoi will dominate 
the cross-scan errors as plotted in Figure IT7I It is for this 
reason that the error in the reference phase in the case of 
Figure DTI has been reduced to 1'. The overall behaviour of 
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Figure 16. This figure shows the mean error in the recovered 
cross-scan position versus the location in phase of the point 
source, for different signal-to-noise ratios. 
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Figure 17. This figure shows the mean error in the recovered 
cross-scan position versus the roll angle, for different signal-to- 
noise ratio point sources. The error on the reference phase in this 
case is 1'. 

the error in the cross-scan position versus roll angle is the 
same in both cases. 

The area available to this method will be the same as 
the area available to the roll angle method as described 
above, since the same criteria on the position of the point 
source relative to the two scans and its detection thresholds, 
as described above, must be met. 

6.3 The number of point sources available to the 
two ring method 

Using the predicted numbers of point sources above the 
given signal-to-noise ratios, as found in Section |1J together 
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Table 5. SNR and the area available for roll angle, p = 0?0. 



SNR 


% area available for p = 0?0 


50 


32.11 


100 


43.02 


200 


53.77 


100 


60.60 



Table 6. SNR versus the number of extragalactic point sources 
available for roll angle, p = 0?0 



SNR 


Total number 
> SNR 


number available for p = 0?0 
> SNR 


50 


59 


18.7 


100 


22 


9.4 


200 


9 


4.8 


400 


5 


3.0 



with the area available to this method found above, it is pos- 
sible to estimate the number of point source detections with 
which the roll angle may be evaluated. The area available, 
for a given roll angle, depends not only on the signal-to-noise 
ratio of the point source, but also on its location in phase as 
seen in Figure |H] It is therefore necessary to know the dis- 
tribution of the point sources in phase before the numbers 
available to this method may be determined. 

6.3.1 Extragalactic point sources 

In the case of extragalactic sources, it may be assumed that 
the point sources are distributed uniformly throughout the 
sphere and hence uniformly in phase. The number of sources 
available to this method therefore, for a given roll angle, 
may be estimated using Figure [T2l The value for percentage 
area accessible to this method, at zero roll angle, is shown 
in Table |SJ using these values the number of point sources 
available to this method may be estimated, and are shown 
in Table El 

The percentage area available decreases with the in- 
creasing magnitude of the roll angle, as seen in Figure 
hence the number of available point sources will also de- 
crease with the increasing roll angle. In the case of zero roll 
angle the number of extragalactic point sources available is 
less than one per week, so it is clear that the evaluation of 
the roll angle will depend predominantly on the galactic and 
not extragalactic point sources. 

6.3.2 Galactic Point Sources 

A uniform distribution in phase is not a good assumption 
in the case of galactic point sources, since their distribution 
in phase will be determined by the position of the galac- 
tic plane in terms of phase. This in turn depends upon the 
scanning strategy used. In order to estimate the numbers of 
galactic point sources available to the method, assumptions 
must be made on the distribution of the point sources in the 
galactic plane and on the scanning strategy employed. The 
numbers of galactic point sources available to the method 
will also depend on the time of year of the observations due 



Table 7. This table displays the number of point sources avail- 
able, for a given roll angle, p, assuming that the spin axis remains 
in the ecliptic plane throughout the mission. The total number 
of point sources available over the course of the mission, assumes 
that the sky is covered twice. This table displays the number of 
point sources available which would provide a measurement of 
the roll angle with an error less than the requirement of 4.4' as 
discussed in Section |21 





p = 0?0 


p = -1?0 


p = +2?0 


Total number 


355.0 


169.1 


288.7 


over mission 








weekly minimum 


1.6 


0.8 


0.7 


weekly maximum 


20.8 


10.8 


20.1 


weekly average 


6.9 


3.3 


5.6 



to constraints on scanning strategy, determined by the rela- 
tive locations of the Earth, Sun, and the spin axis position. 

In order to investigate these dependencies, the phase 
coverage of the galactic plane was evaluated on a weekly 
basis for a given scanning strategy. The number of point 
sources available to this method may then be found, both 
on a weekly basis and over the course of the whole mission. It 
was found that the distribution assumed in galactic latitude 
had little effect on the number of point sources available, 
whereas the assumed distribution in galactic longitude had 
a large effect on the variation of the number of point sources 
available per week. 

In order to obtain the most reliable estimates of the 
number of point sources available on a weekly basis; the 
galactic point sources found in Section |1J were binned into 
10° bins in galactic longitude, and within these bins they 
were assumed to be uniformly distributed between -5° and 
+5° in galactic latitude. 

The first scanning strategy considered, was the base- 
line strategy of the spin axis remaining in the ecliptic plane 
throughout the course of the mission, with the nominal 1° 
motion in ecliptic longitude per day, maintaining the spin 
axis in the sunwards direction. Here the boresight angle is 
defined as the angle between the negative, anti-sunwards di- 
rection, spin axis and the line of sight. Table shows the 
number of point source available, for different values of the 
roll angle, from which the roll angle may be successfully eval- 
uated with an error below the required limit, as discussed in 
Section|5| The table shows the minimum, maximum and av- 
erage number of sources available in the course of one week 
over the course of the mission together with the total num- 
ber of point sources available over the whole mission; here 
the mission is assumed to last for the nominal one year with 
two sky coverages. Table Q also shows the effect of the de- 
crease in the area available to this method with increasing 
magnitude of the roll angle, as seen in Figure FHH with the 
decrease in the number of available point sources for the roll 
angle values of p — — 1?0 and p = +2?0. 

Sinusoidal scanning strategies may also be considered. 
These scanning strategies may be defined by the path of the 
spin axis position (/3 Z ,A Z ), over the course of the mission, 
given by: 

A z = A, 

I3 Z = Asm{nX + <j)), (25) 
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Table 8. The variation of the weekly minimum value over various 
sinusoidal scanning strategics, for different roll angles. 



Roll angle 


Weekly minimum 




(°) 


minimum maximum 


average 


-1.0 


0.4 1.4 


0.8 


0.0 


1.1 2.1 


1.6 


+2.0 


0.4 1.8 


0.8 



Table 9. The variation of the average weekly value over various 
sinusoidal scanning strategies, for different roll angles. 



Roll angle 


Weekly average 




(°) 


minimum maximum 


average 


-1.0 


2.9 4.4 


3.5 


0.0 


6.2 7.7 


6.9 


+2.0 


4.9 6.6 


5.6 



where A is the amplitude of the motion out of the ecliptic 
plane, A=0° is equivalent to the above scanning strategy 
with the spin axis remaining in the ecliptic plane throughout 
the whole mission. <f) defines where the spin axis crosses the 
ecliptic plane and n is the number of periods of the sinusoidal 
motion over the course of the mission. 

The parameters, A, n and <f> were allowed to vary in 
order to investigate possible scanning strategies which will 
maximise the number of point sources available, in partic- 
ular those scanning strategies with larger weekly minimum 
values. The value of A is constrained to be ^ 10° as the spin 
axis must remain within 10° of the Sun, this constraint is 
due to the physical size of the solar panels and the need to 
minimise stray-light. 

It was found that the value of (j> had the most effect on 
the number of point sources available. The values of A and n 
also have an effect on the number of point sources available, 
but the values which maximise the number available depend 
on the value of the roll angle. Given that the value of <j> 
cannot be predetermined as it too is constrained, by the 
Earth aspect angle which must remain ^ 15° within the L2 
orbit, it is not possible to recommend a particular scanning 
strategy for the benefit of the roll angle reconstruction. 

Tables |H| and [5] show the variation in the numbers of 
point sources available over the range of sinusoidal scanning 
strategies investigated. Table |H] shows the variation in the 
value of the weekly minimum for different values of roll an- 
gle, from this table it may be seen that as the magnitude 
of the roll angle increases so does the likelihood of there be- 
ing an occasional week in which there will not be a suitable 
point source available to extract the current value of the roll 
angle. 

6.4 The impact of the error in the recovery of the 
effective boresight 

The discussion so far has focused on the extraction of the 
roll angle, assuming that it is allowed to consume the total 
error budget for the pointing reconstruction. Considerations 
of the error in the reconstruction of the effective boresight, 
a, must now be included in the error budget. The error in 



Table 10. This table shows the requirements on the error in the 
roll angle, given an error in the effective boresight angle. The 
error in the effective boresight assumes an error in the cross-scan 
position of the point source of <r c =3", an error in the spin axis 
position of 1" and an error in the point source position as given 
below. The requirements on the error in the roll angle are given 
for the two cases of 70% and 100% sky coverage. 
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J sky 


0.7 f sky = 1.0 




(") 


(') 


(') 


(") 
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4.6 


4.0 


3.6 


5 


4.4 


3.7 


5.6 


6 


4.1 


3.4 


7.3 


7 


3.7 


3.0 


8.9 



Table 11. The variation of the weekly minimum value over var- 
ious sinusoidal scanning strategies, for different requirements on 



the roll angle 


error, for p 


= 0?0. 




Roll angle 


Weekly minimum 




Error ^ 


minimum 


maximum 


average 


4.4' 


1.1 


2.1 


1.6 


4.0' 


1.0 


2.2 


1.5 


3.7' 


1.0 


1.8 


1.3 


3.4' 


0.7 


1.7 


1.1 


3.0' 


0.6 


1.5 


0.9 



the effective boresight depends on the error in the angle 
found between the spin axis position and a point source 
of known position and the error in the cross-scan position 
found for the point source. Only the error in positions in the 
cross-scan direction will have an effect on the error in the 
effective boresight. If it is assumed that the positional errors 
are the same in both the scan and cross-scan directions then 
an expression for the error in the effective boresight, a a , may 
be defined as: 

fa = g (°Pt + a sa) + °c, (26) 

where cr^ is the error in the cross-scan position of the point 
source and a p t and a sa are the radial errors in the position 
of the point source and the spin axis, respectively. 

Table fTTTI shows the effect of including the errors in the 
effective boresight on the requirements for the error in the 
roll angle in order that their combined effect on the point- 
ing error is less than the requirements found in Section |5] 
Included in this table are the requirements on the radial 
error in the point source position if it is assumed that the 
error in the cross-scan position found for the point source is 
3", which is the worst case value for a point source with a 
signal-to-noise ratio of 50 as seen in Figure [Trjl 

Tables 1111 and 1121 show the effect of the increased re- 
quirements on the accuracy of the roll angle, on the vari- 
ations in the weekly minimum and average values of the 
number of available point sources over the various sinusoidal 
scanning strategies. In Table 1111 the minimum value of the 
weekly minimum drops below one for the requirement that 
the error in the roll angle, a p , is ^ 3.4. The average value 
of the weekly minimum is less than one for the requirement 
that Op is ^ 3.0. In this table the value of the roll angle is 
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Table 12. The variation of the weekly average value over various 
sinusoidal scanning strategies, for different requirements on the 
roll angle error, for p = 0?0. 



Roll angle Weekly average 

Error ^ minimum maximum average 



4.4' 


6.2 


7.7 


6.9 


4.0' 


5.7 


7.4 


6.4 


3.7' 


5.3 


6.8 


5.9 


3.4' 


5.0 


6.2 


5.5 


3.0' 


4.4 


5.7 


4.9 



zero, which we have earlier established is the value of the 
roll angle with the greatest number of point sources avail- 
able for its recovery. It is therefore desirable that the limit 
of the error on the roll angle is as large as possible in order 
to increase the range in the values of roll angles for which 
there may be a weekly reconstruction of the roll angle. 



7 DISCUSSION 

We have shown here, in Section 13.1.11 that the roll angle 
cannot be evaluated using data from a single ring. The cross- 
scan position of a point source cannot be evaluated before 
the roll angle has been found as the cross-scan positions 
of the detectors are unknown. The cross-scan position of a 
point source can, however, be evaluated using data from a 
single ring, once the roll angle is known. The cross-scan posi- 
tion of a point source may also be evaluated by an extension 
to the method used to reconstruct the roll angle from two 
neighbouring rings, as discussed in Section [3.1.21 The eval- 
uation of the cross-scan position using two rings is immune 
to the relative calibration errors between the detectors, but 
the area available to this method is the same as the area 
available for the reconstruction of the roll angle; where the 
area available depends on both the signal-to-noise ratio of 
the point source and on the actual value of the roll angle. In 
comparison the area available to the extraction of the cross- 
scan position from a single ring is only dependent on the 
value of the roll angle, through its effect on the cross-scan 
positions of the detectors. In order to extract its cross-scan 
position, the point source must lie within the range of the 
cross-scan positions of the detectors. This criterion leads to 
an opposite dependence on the roll angle, as compared to 
the two ring method; the area available increases with the 
magnitude of the roll angle, with the minimum value occur- 
ing at zero roll angle. 

The main components of the error in the effective bore- 
sight are the error in the recovered cross-scan position and 
the error in the position of the point source. The error in the 
effective boresight affects the requirement on the roll angle 
error in order that the resultant error in the pointing meets 
the requirements as found in Section|5| The less stringent the 
requirements on the roll angle error, the greater the range 
in the roll angle which may be reliably reconstructed on a 
weekly basis. It is therefore possible to place requirements 
on the radial error in the position of the point source used in 
the evaluation of the effective boresight. Using typical values 
of the error in the recovered cross-scan position, it is found 
that in order to meet the pointing requirement, the radial 



position error should be less than 10". Ideally the radial er- 
ror should be ~3" in order to reduce the requirements on 
the recovery of the roll angle. 



8 CONCLUSIONS 

We have investigated methods for the evaluation of param- 
eters required for the pointing reconstruction of the Planck 
satellite. The effect of pointing errors on the recovered Ces 
has been assessed and requirements on the accuracy of the 
pointing reconstruction found, such that the effect of ran- 
dom pointing errors is less than the unsubtractable noise. 
The accuracy at which the roll angle is required is deter- 
mined so that the uncertainty in the positions of the detec- 
tors due to the error in the roll angle is less than the limit 
on the pointing error. 

It is found that it is impractical to determine the roll 
angle from data from a single ring, as this would require an 
extremely high signal-to-noise ratio point source to success- 
fully determine the roll angle. An alternative method using 
two neighbouring rings was developed, which allows the suc- 
cessful extraction of the roll angle to the required accuracy 
from a single point source, with a signal-to-noise ratio ~ 50. 
This method also has the additional advantage of using the 
ratio of the amplitudes of a point source as seen by the same 
detector on different rings, hence this method does not re- 
quire the relative calibration of the detectors. This method, 
however, does rely upon the assumption that the roll angle 
is a slowly varying parameter, and it should be remembered 
that any abrupt changes to the inertia tensor will invalidate 
this assumption. Once the roll angle is determined the cross 
scan position of any point source may be assessed. The effec- 
tive boresight angle may then be evaluated using the cross 
scan position of a point source of known position. 

Due to the depletion of consumables over the course of 
the mission and the resultant variation in the values of the 
pointing parameters, these parameters should be evaluated 
on a regular basis. Whether a method may be applied suc- 
cessfully on a regular basis, will depend on the availability 
of suitable point sources over the course of the mission. In 
order to investigate this availability, the number of galac- 
tic and extragalactic point sources detectable by the Planck 
mission was determined. It was found that there is an in- 
sufficient number of suitable extragalactic point sources to 
evaluate these parameters, the roll and effective boresight 
angles, on a weekly basis. Hence the number of galactic 
point sources available is crucial. The total number of galac- 
tic point sources observable by Planck over the course of the 
mission, and more importantly the variation in the number 
observable on a weekly basis, depends upon the scanning 
strategy used. However, it was found that it was not possi- 
ble to recommend a scanning strategy which will maximise 
the number of galactic point sources available. This is due 
to other constraints on the scanning strategy such as main- 
taining the Earth aspect angle ^15° within the L2 orbit. 

The number of point sources available is also found to 
depend on the value of the roll angle. The method, pre- 
sented here for the reconstruction of the roll angle from two 
neighbouring rings, is best suited for the evaluation of small 
magnitude roll angles. The larger the roll angle, the smaller 
the number of point sources available and the greater the 
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chance of weeks with no suitable point sources. For roll an- 
gles ~ +2° and ~ — 1° it is likely that there will be an 
occasional week in which there is no suitable point source 
to evaluate the roll angle. As the magnitude of the roll an- 
gle increases above these values the number of suitable point 
sources diminishes until ~ +4° and ~ —3° where the evalu- 
ation of the roll angle using the method presented here is no 
longer possible. This places a manufacturing requirement on 
the alignment of the focal plane with respect to the nominal 
scan direction. 

While the evaluation of the roll angle does not require a 
known point source, the reconstruction of the effective bore- 
sight requires the position of the point source to be known. 
A known point source catalogue will therefore be required 
for the evaluation of the effective boresight angles, with po- 
sitional errors less than 7" and ideally less than 2", in order 
to meet the accuracy requirements on the pointing parame- 
ters. 

Given the success of this method depends on the num- 
bers of bright point sources available, and for the reconstruc- 
tion of the effective boresight on their accurate positions, it 
is important that the IRAS sources which are found here 
to be potentially suitable for the reconstruction of these pa- 
rameters are observed prior to launch, to confirm their point 
source nature, to improve upon their flux estimates and ac- 
quire accurate positions. Any point source ^ 50 times the 
noise expected in a binned ring in the 857 GHz channel 
is potentially useful in evaluating these parameters. Given 
the expected noise level this corresponds to point sources 
> 14.2 Jy at 857 GHz. 

The method presented here will allow the weekly mon- 
itoring of the roll angle, given the number of point sources 
determined here, under the expectation of small roll angles 
and a slowly varying inertia tensor. When combined with a 
catalogue of accurate point source positions, this method 
will also allow the weekly determination of the effective 
boresight angle. It is expected that the determinations of 
these two angles from the science data will be used to de- 
termine their systematic offsets from the equivalent angles 
determined more accurately by the Star Tracker. 
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